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Based on a semi-empirical generalized Anderson-Newns model we construct a pseudo-particle 
description for electron emission due to de-excitation of metastable molecules at surfaces. The 
pseudo-particle approach allows us to treat resonant charge-transfer and Auger processes on an 
equal footing, as it is necessary when both channels are open. This is for instance the case when 
a metastable N2( 3 E^) molecule hits a diamond surface. Using non-equilibrium Green functions 
and physically motivated approximations to the self-energies of the Dyson equations we derive a 
system of rate equations for the probabilities with which the metastable N2( 3 E„) molecule, the 
molecular ground state N2( 1 S+), and the negative ion Nj( 2 n g ) can be found in the course of the 
scattering event. From the rate equations we also obtain the spectrum of the emitted electron and 
the secondary electron emission coefficient. Our numerical results indicate the resonant tunneling 
process undermining the source of the Auger channel which therefore contributes only a few percent 
to the secondary electron emission. 

PACS numbers: 34.35. +a, 34.70.+e, 79.20.Fv, 79.20.Hx 



I. INTRODUCTION 

Charge exchange processes during atom-surface or 
molecule-surface collisions have been the subject of in- 
tense scientific research during the last decades-^. This 
type of surface reactions is of fundamental interest. It 
represents a quantum-impurity problem where a finite 
many-body system with discrete quantum states couples 
to an extended system with a continuum of states which 
essentially acts as a reservoir for electrons. Under ap- 
propriate conditional^ such an arrangement gives for in- 
stance rise to the Kondo effect^ originally found in metals 
containing magnetic impurities or to Coulomb blockades 
as it is discussed in nanostructurea^l 

Besides of being a particular realization of a quantum 
impurity problem, atom/molecule-surface collisions are 
also of technological interest, especially in the field of 
bounded low-temperature plasmas, where this type of 
surface collisions is the main supplier of secondary elec- 
trons which in turn strongly affect the overall charge 
balance of the discharge-. In dielectric barrier dis- 
charges, for instance, secondary electron emission de- 
termines whether the disc harge operates in a filamen- 
tary or a diffuse mode^^. Only the latter mode is use- 
ful for surface modification. Controlling the yield with 
which secondary electrons are produced is thus of great 
practical interest. This applies even more so to micro- 
discharges^ where the continuing miniaturization gives 
charge-transferring surface reactions more and more in- 
fluence on the properties of the discharge. 

Depending on the projectile and the target, secondary 
electron emission usually occurs either in the course of 
a resonant tunneling process or an Auger transition. In 
some situations, however, both transitions may be en- 
ergetically allowed and hence contribute to the yield 
with which electrons are released. The interplay of the 
two reaction channels has therefore b een studied in th e 
pastP^H^. Starting with the work of Alvarez et al.'^ a 



detailed theoretical analysis of the interference of Auger 
and resonant tunn eling pro cesses has been given by Gold- 
berg and coworker j 18 l 2 ° l 21 l. Their results for H + and He + 
indicate the Auger channel to be active only close to 
the surface whereas the resonant channel is already ef- 
ficient at rather large projectile-surface distances. When 
both channels are coupled together the dynamics of the 
system is hence controlled by the resonant channel as 
it destroys the initial species before the Auger channel 
can become operative. The Auger channel is therefore 
strongly suppressed in the coupled system albeit the in- 
dividual efficiencies of the reaction channels are compara- 
ble. Onufriev and MarstorP^ also investigated the inter- 
play of tunneling and Auger processes for the particular 
case of Li(2p) atoms de-exciting on a metallic surface. 
Using a sophisticated many-body theoretical description 
of the scattering process they concluded that depending 
on the model parameters the two de-excitation channels 
interfere either constructively or destructively. 

Whereas the previous studies focused on atomic pro- 
jectiles we will investigate in the following the interplay 
of Auger and resonant tunneling processes for a molecu- 
lar projectile. More precisely, we will analyze how these 
two processes affect secondary electron emission due to 
de-excitation of metastable molecules. Neutralization of 
molecular iomP2H23] w iH not be discussed. A particu- 
larly interesting case is the de-excitation of metastable 
N2( 3 S^) becaus e this molecule de-excites in two primary 
reaction channels^ 25 * 26 !. On the one hand, there is the two- 
step resonant charge transfer (RCT) reaction 



eg+N 2 ( d £+) 



4N 2 -( 2 n ff ) 



^N 2 (^+) + e f , (1) 



where e- and e- denote an electron within the surface 
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and a free electron, respectively. In this process the 
metastable N2( 3 £+) molecule first resonantly captures 
an electron from the surface to form the intermediate 
negative ion shape resonance ( 2 II g ) which then decays 
into the ground state N2( 1 Si") by resonantly emitting an 
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electron. The decay of the negative ion can be either due 
to the ion's natural life time or due to the ion surface 
interaction. In addition to the RCT channel, there exists 
an Auger de-excitation reaction also known as Penning 
de-excitation 



eg+N 2 ( J £+: 



+ N 2 (^ 



(2) 



Here the molecule non-resonantly captures a surface elec- 
tron and simultaneously releases another electron. De- 
pending on the surface band structure both processes 
may be possible at the same time. This is for instance 
the case for diamond as it possesses a rather wide valence 
band. 

In our previous worlP^U we investigated the two reac- 
tion channels separately by a quantum-kinetic approach 
and a rate equation technique. The molecule was in 
both cases treated as a semi-empirical two-level system 
corresponding to the 2tt u and 2-7r g molecular orbitals, 
which are the two molecular orbitals whose occupancies 
change during the de-excitation process. Spectator elec- 
trons not involved in the processes were neglected. De- 
pending on the process the two levels denoted the up- 
per and lower ionization levels of either the metastable 
molecule (Auger de-excitation) or the negative ion (reso- 
nant charge-transfer) . The advantage of a semi-empirical 
model is that it is based on a few parameters which are 
relatively easy accessible, either experimentally or theo- 
retically. The difference of the two parameter sets, which 
arises from intra-molecular Coulomb correlations not in- 
cluded in the model, is not a problem as long as the two 
processes are treated separately. A simultaneous treat- 
ment of them requires however a way to implement both 
parameterizations within a single Hamiltonian. 

A way to overcome the problem would be of course to 
set up a more general projectile Hamiltonian, including 
active and spectator electrons and their mutual Coulomb 
interactions. Treating these intra-molecular Coulomb 
correlations explicitly is however extremely demanding. 
It embraces a full quantum-chemical description of the 
approaching molecule which cannot easily be adapted 
from one projectile to another. Since we are primar- 
ily interested in developing models and tools to be used 
for the description of secondary electron emission from 
plasma walls easy adaptation from one projectile-target 
combination to another is however an important crite- 
rion for us. We stay therefore within the limits of a semi- 
empirical two-level system and use instead projection op- 
erators and two auxiliary bosons to assign and control 
the level energies. With these constructs it is possible 
to formulate a Hamiltonian containing both channels - 
and (2j| - without introducing intra-molecular interac- 
tions. The projection operators allow us to assign differ- 
ent parameterizations to the two levels, depending on the 
occupancy of the system, whereas the auxiliary bosons 
enable us to mimic the intra-molecular Coulomb corre- 
lations which need to kick in to make the two tunneling 
processes involved in ([I]) resonant. Expressing the projec- 
tion operators in terms of pseudo-particle operators with 



boson or fermion statistics opens then the door for em- 
ploying non-equilibrium Green functions to derive from 
the Hamiltonian quantum-kinetic equations for the prob- 
abilities with which the molecular states can be found in 
the course of the collision. 

The strength and flexibility of the pseudo-particle or 
slave field approach, originally developed by Coleman^ 
in the context of the infinite-?/ Anderson Hamiltonian, 
has been demonstrated many ti mes for Anderson-type 
and Anderson-Newns-type models^^HH] yy e apply this 
method to a generalized Anderson-Newns model describ- 
ing the coupling of different molecular configurations to 
a solid. It is of the type but not ide ntical to the model in- 
troduced by Marston and coworker d 14 * 34 l The derivation 
of the quantum-kinetic equations for the pseudo-particle 
propagators with the subsequent reduction to the rate 
equations for the occupancies of the molecular pseudo- 
particle states follows the work of Langreth and cowork- 
er j5 | 3i | a s a result we obtain rate equations describ- 
ing the de-excitation of N 2 ( 3 E+) in situations where the 
RCT and the Auger channel are simultaneously open. 
In the absence of the Auger channel the rate equations 
reduce to the ones we derived intuitively before for the 
isolated RCT channeP^. Applying the model to the par- 
ticular case of a diamond surface shows the resonant 
process dominating the Auger process. The overall sec- 
ondary electron emission coefficient due to de-excitation 
of N 2 ( 3 £+) at a diamond surface is on the order of 10 _1 . 

The outline of the rest of the paper is as follows. In 
Sec. [IT] we describe the semi-empirical model on which 
our investigation of the de-excitation process is based. 



Thereafter, we explain in Sec. Ill the pseudo-particle rep- 
resentation. Afterwards, we conduct in Sec. |IV| a second 
order quantum-kinetic calculation on top of the pseudo- 
particle model. In Sec. [V] we introduce a physically mo- 
tivated semi-classical approximation that allows us to re- 
duce the set of Dyson equations to a set of rate equations. 
Finally, we present in Sec. [VI] the results for the diamond 
surface and conclude in Sec. VII| with a brief summary 
of the main points of the work. Appendix [A] lists the 
Langreth- Wilkins rules^ as used in our calculation and 
Appendix [B] collects the second order Dyson equations 
for the molecular Green functions. 



II. MODEL 

The interacting molecule-surface system is character- 
ized by three different types of electronic states: bound 
and unbound molecular states and states withi n the 
solid surface. In the spirit of our previous worlPZM] 
we restrict the attention to those states whose occupan- 
cies change during the molecule-surface collision. For 
these states and the coupling between them we con- 
struct a semi-empirical model. Its matrix elements can 
be either obtained from quantum-mechanical calcula- 
tions based on particular assumptions about the electron 
wave functions and/or experimentally measured ioniza- 
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N +( 2 n u ) N 2 ( 1 Sg ) n 2 ( 3 ^) N^( a n tf ) 

FIG. 1. Correspondence between the occupation of the two 
level system and the molecular states. Here eo and ej denote 
the energies of the levels and 1, respectively. 



tion energies, electron affinities, surface response func- 
tions, and electron tunneling rates. Since we are pri- 
marily interested in the quantum-kinetic handling of the 
semi-empirical model we pursue for simplicity the former 
route. A more realistic parameterization of the model is 
however in principle possible. 

We treat the relevant bound states of the nitrogen 
molecule in terms of a two-level system consisting of a 
ground state level "0" and an excited level " 1" . Within a 
linear combination of atomic orbital (LCAO) description 
of the molecule these two levels represent the nitrogen 
molecule's 2tt u and 2ir g orbitals. Each of the two or- 
bitals can carry four electrons. We neglect however the 
three electrons in the 2ir u orbital and the three holes in 
the 2ir g orbital which are not directly involved in the 
de-excitation process. They act only as frozen-in spec- 
tators. For the same reason we neglect the electron spin 
and treat the magnetic quantum number to = ±1 as an 
initial parameter. Hence, both levels of our model carry 
at most one electron. 

The two-level system represents any of the molecular 
states depicted in Fig. [I] The positive ion N^ 2 !^), the 
ground state N 2 ( 1 Sg ), the metastable state N 2 ( 3 £+), 
and the negative ion N^"( 2 n s ). Depending on the par- 
ticular occupation the energies eo and £i correspond 
therefore to the ionization energies of different molecu- 
lar states. Due to intra-molecular Coulomb in teract ions 
these ionization energies are in general different 36 ! 37 !. We 
have to allow therefore eo and e\ to depend on the occu- 
pancy of the two-level system, that is, on the molecular 
state it is supposed to represent. In addition, the ion- 
ization energies are also subject to the surface's image 
potential Vi(z) and thus vary with time as the molecule 
moves with respect to the surface. Using the analysis 
presented in Ref . [55] we find 



£o g {z) =s x g -V i (z) , 
eu(z)=Eg-Vi{z) , 
eo-(*)=e~ +V5(*)> 



£l-{z) 



Vi(z) , 



(3a) 

(3b) 
(3c) 
(3d) 



where the subscripts g, * and — signal the dependence 
of the energy levels on the molecular state denoting, 
respectively, the ground state molecule, the metastable 
molecule, and the negative ion. The unperturbed molec- 
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vacuum level 



Ae v 



E 



(c) 



Sl-/q{z) 



£o-0) 



Diamond 



FIG. 2. Energy scheme for the case of a metastable N2( 3 E+) 
molecule interaction with a diamond surface. The Auger 
de-excitation channel (a) is depicted on the left hand side 
whereas the RCT electron capture (b) and the subsequent 
RCT electron release (c) are shown on the right hand side. 
The latter transition is drawn with a dashed line because the 
electron emission actually occurs at the same z position and 
hence energy. It is resonant. Therefore the emitted electron 
will only reach the shown location over time. The drawing is 
to scale in terms of energy units. 



ular energies are given by 

e- = -17.25 eV 
e 2° = -14.49 eV 



..36 37 



e^ = -9.67eV 
eT = 1.18 eV . 



(4) 



The overall energy scheme of the coupled molecule- 
surface system is sketched in Fig. [2] for the particular 
case of a diamond surface. As can be seen the positive 
ion is neither involved in the RCT nor the Auger process. 

For the image potential we employ for simplicity the 
classical expression, 



Vt(z) « - 



1 



e b r + 1 167reo z 



(5) 



with e b r standing for the surface's static bulk dielectric 
constant. Close to the surface the image potential is how- 
ever truncated according to Vi(z c ) = Vq where Vq is the 
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depth of the potential barrier confining the electrons of 
the solid participating in th e de-e xcitation process. As 
in our previous investigations^!^ W e describe the solid 
by a step potential of depth Vq. For a metallic surface^ 
the step depth is the width of the conduction band Aec, 
that is, the sum of the work function $vv and the Fermi 
energy e p whereas for a dielectric surface^ it is the sum 
of the width of the valence band Aey, the energy gap 
e g , and the electron affinity e a which can be positive or 
negative. 

In both the RCT and the Auger channel the emitted 
electron stems from the molecule and, thus, the emission 
proceeds into the molecular continuum states. We model 
the latter as free electron states moving along with the 
molecule and label them with q. Electrons residing in 
those states are also affected by the image potential and 
their energy is thus given by 



e q ~(z)=ef+V i (z) = ^ r +V i (z). (6) 



The wave function of the emitted electron is a two-center 
Coulomb wave for the Auger channeP^ and a plane wave 
for the RCT channel 2 -. The latter is a special case of 
the former which holds for zero effective nucleus charge. 
Since in the RCT channel the emitted electron leaves 
a neutral molecule behind the plane wave is the suitable 
choice for this channel. In the Auger channel however the 
emitted electron feels the residual two-center Coulomb 
attraction of the ion core. The two-center Coulomb wave 
takes this effect into account. To complete the descrip- 
tion of the model we note that we use LCAO molecular 
wave functions for the two-level system and the eigen- 
states of the step potential for the wave function of an 
electron form the solid. Explicit expressions for the wave 
functions and the details of the calculation of the Auger 
and tunneling matrix elements are given in Refs. 1271 and 
El 

We now cast the semi-empirical model just described 
into a mathematical form. Introducing projection oper- 
ators 



P" " 1 = |noni)(noni| (7) 



projecting onto states of the two-level system with no = 
0, 1 electrons in the lower and ri\ = 0, 1 electrons in the 
upper state, the transitions shown in Fig. [2] give rise to 



a generalized Anderson-Newns modeP^, 
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+ J2 pnoni (er il {t)clc + eT ni {t)c\ Cl ) 

+ Y,([P 01 + P 11 ]V,(t)cUlc + h.c) (8) 

k 

+ J2([P 1Q + P 11 ]v q ~(t)cib{c 1 + h.c) 

kq 

where ej° = e 0g , E® 1 = £u, e^ 1 = e -, and e} 1 = £i_. 
The remaining energy levels need not be specified fur- 
ther. They drop out in the course of the pseudo-particle 
representation presented in the next section. The two 
auxiliary bosons 6q and b± mimic the intra-molecular 
Coulomb correlations which make the two steps of the 
RCT channel (IT]) resonant. This can be accomplished by 
setting ujq = — e^ 1 and cjj = ej 1 — e^ . The initial 
energy of the tunneling electron is then resonant respec- 
tively with the lower and the upper level of the negative 
ion. The rest of the Hamiltonian is wr itten in the nota- 
tion we used in our previous worfP^U. 

The time dependence of the Hamiltonian arises from 
the trajectory of the molecule's center of mass R(t). For 
simplicity we assume normal incidence described by the 
trajectory R(t) — (vq \t\ + zo)e z , where vq is a constant 
velocity and z is the molecule's turning point which 
can be calculated from a Morse type interaction poten- 
tial. Employing the trajectory the z dependence of 
the molecular energies ^ and the energy of the emitted 
electron (I6J) transforms into a time dependence. Simi- 
larly, the^uger matrix element Vg- and the two reso- 
nant tunneling matrix elements Vg and Vq acquire also a 
time-dependence. 

The quantum-kinetic calculation presented below 
treats the matrix elements of the Hamiltonian as param- 
eters. We are thus not restricted to the specific approx- 
imations to the matrix elements derived in our previous 
worlpZISS]. We could as well use matrix elements from ab- 
initio calculations or experimental measurements. The 
natural decay of the negative ion, described by the rate 
r„ = 1/t„ with t„ the natural lifetime of the negative 
ion, is not included in the Hamiltonian ([8]). It will be 
inserted at the end into the final set of rate equations for 
the molecular occupancies. 

III. PSEUDO-PARTICLE REPRESENTATION 

The projection operators ([7| permit us to describe the 
transitions (fTl) and pi) by a single Hamiltonian. Dc- 
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pending on the process and thus the occupancy of the 
molecular levels different matrix elements can be assigned 
to the Hamiltonian. The projectors also guarantee that 
the N^( 2 n u ) state never occurs. In other words, they 
ensure that the occupancies of the two molecular levels 
never vanishes simultaneously. For instance, an electron 
residing in the upper level can only be released when an 
electron has been captured in the lower level. 

A drawback of the projection operators is that they are 
not suitable for a diagrammatic treatment which on the 
other hand is a powerful tool to set up quantum-kinetic 
equations. To remedy this drawback we now employ a 
pseudo-particle approach to express the Hamiltonian Q 
in terms of slave nelda^sHU]. The starting point for this 
procedure is the completeness condition, 



where we have introduced the usual pseudo-particle num- 
ber operator Q to encapsulate the constraint. 

Formally, the auxiliary fcrmion and boson operators 
are pseudo-particle operators creating and annihilating 
molecular configurations. The constraint ensures that at 
any time only one of the four possible molecular con- 
figurations is present in the system. The occupancy 
of a molecular pseudo-particle state is thus at most 
unity. Hence, it represents the probability with which 
the molecular configuration it describes appears in the 
course of the scattering event. 



Inserting the decomposition (12 1 into the Hamiltonian 
([8]), making the identifications (13), and collecting only 
terms which are in accordance with ( 14 ) is straightfor- 
ward. The result is 



|00)(00| + |10)(10| + |01)(01| + |11)(11| 



1 



(9) 



which expresses the fact that the molecule can be only 
in either one of the configurations depicted in Fig. [I] 



Introducing pseudo-particle operators c+, c g , cl, and c 
that create the positive ion, the ground state molecule, 
the metastable molecule, and the negative ion from an 
abstract vacuum state, 



c^vac), 1 10) 



|00) 

|01) = 4 |vac), |ll)=ct|va 
the completeness condition ^ becomes 



(10a) 
(10b) 



+ C' g C g + C\ C„ + CL C_ = l 



Using ([9]) and ( 10 ) the operators creating and de- 
stroying an electron in the two states of the two-level 
system can then be written as 



c 


= c * 1 = 


|00><10| - 


|oi><ii 


c f 
<-o 


= c * 1 = 


|10><00| - 


|ii><oi 


c l 


= c x * 1 = 


|00)(01| + 


lioxii 


°1 


= c\ * 1 = 


|01)(00| + 


liiXio 



(12b) 



\ = c\c, +cjc_ , (12c) 



(12d) 



cl c , 



cl c„ 



In order to satisfy the anti-commutation rela- 
cj^ one has to define c |11) = — 1 01) 

and Co|01) = -|H) (see also Ref. 138]). The anti- 
commutator relations then reproduce the completeness 
condition 



the 



111 when cither the c s /* are bosonic and 



7+ 



are fcrmionic or the 



z g/* 



are fermionic and 



the C.-/+ are bosonic. Without loss of generality we 
choose c g and c* to be bosonic and declare the labeling 
conventions 



„(t) 

„(t) 



f- 



(t) 



„(t) 



L(t) 



(f) 



The constraint (11) reduces then to 

bt K + /!/_ + /I A 



Q 



b l b 9 



1 , 



(13a) 
(13b) 

(14) 



+ e g (t)blb g + e m (t)btb m +e_(t)f*_f_ 



Y J {v i {t)c\blb\f_+h.c. 

k 

Y,(y^t)c\b\b\f_ + h.c 

J2( v ^ c A b l b * + h - c ) . 



(15) 



kq 



(11) where we introduced the abbreviations 



£ 9 =£o°: £*=4 1 , e_=el 1 +e\ 1 . (16) 

Notice, no term in the Hamiltonian contains the opera- 
tor f + or its adjoint. This is by construction because the 
positive ion N^( 2 n u ) is not involved in the two transi- 
tions the Hamiltonian is supposed to model. The phys- 
ical meaning of the various terms in the Hamiltonian is 
particularly transparent. Consider for instance the last 
term describing the Auger de-excitation. A metastable 
molecule and an electron from the surface disappear 
while a ground state molecule and an Auger electron are 
created. 

The operators /_ and b^,/ g comply to standard commu- 
tation and anti-commutation relations. It is thus possible 
to conduct a non-equilibrium diagrammatic expansion 
of the interaction terms in (15). The Hamiltonian (15) 



conserves the pseudo-particle number Q. The quantum- 
kinetic equations however may contain terms which vi- 
olate the constraint. The projection onto the physical 
subspace with Q = 1 needs to be therefore carried out 
explicitly. For this purpose we employ in the next sec- 
tion the Langreth-Nordlander projection technique^. 



IV. QUANTUM-KINETICS 

To start the quantum-kinetic calculation we first de- 
fine a contour-ordered fermion Green function G_ for 
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the negative ion and contour-ordered boson Green func- 
tions -B*/ s /o/i for the metastable molecule, the molecu- 
lar ground state, and the t wo auxiliary bosons, respec- 
tively. Using the notation of Langreth and Nordlandei' ^ 
we write 



^G_(^,^')-(T c /_(^)/!(^')), 

iB l (t,t') = (T c b l (t)bl(t')) , 



(17a) 
(17b) 



where I = *,<?, 0, 1, and define the analytic pieces Gl 

> 

and by 



iG-{t,t') = e c (t-t')G>(t,t') 

-e e (t'-t) G<(t,t') 
iB l (t,t') = e c (t-t')B>(t,t') 

+ e c (t'-t) B<(t,t') 



(18a) 
(18b) 



The time-ordering operator Tq and the 0c function are 
defined on a complex time contour. The associated re- 
tarded Green functions G R and BP are given by 

iG R {t,t') =Q(t-t')(G>(t,t')+G<(t,t')) , (19a) 
iB R (t,t') = G(t-t')(B>(t,t')-B<(t,t')) , (19b) 



k+q 



B*+g(t 1 ,t 2 ) 




(* a ,ti) 






v iqM 



whereas the advanced functions G^ 4 and i3, may be con- 



structed from ( 19 1 using the relations 



G*{t,t') = [G R {t',t)] 
Bf(t,t')=[B R (t',t)]' 



(20a) 
(20b) 



In contrast to our previous worlP^U, we no longer use 
Keldysh's matrix notation. Langreth's notation is more 
convenient. It enables us to use the powerful Langreth- 
Wilkins rules^ which lead more directly to the quantum- 
kinetic equations for the molecular occupancies, that is, 
the probabilities with which the molecular configurations 
involved in the de-excitation process appear in the course 
of the scattering event. 

In order to calculate the self-energies we truncate the 
diagrammatic expansion beyond the second order and 
employ the self- consistent non-crossing approximation^. 
The diagrams for the fermionic self-energy £_ and the 
bosonic self-energies II* and IT ff are shown in Fig. [3] 
Mathematically they read 



+ ia q -(t 1 ,t 2 ) B g (t 1) t 2 ) , 

n*(ii,i 2 ) = -to-jj(t2,ti)G_(ti,t 2 ) 
+ a^(t u t 2 ) B g (t u t 2 ) , 

n g (h,t 2 ) = -^ ( 7 ^ (^ 2 ,^l)G_(^l,^ 2 ) 



(21a) 
(21b) 
(21c) 



FIG. 3. Second order self-energy diagrams in self-consistent 
non-crossing approximation. Depicted are from top to bot- 
tom the self-energy of the fermionic level — iE_(ti,i2), the 
RCT component of the self-energies of the bosonic lev- 
els — ill*(ti,ta) and — iH g (ti,t2), the Auger component of 
the metastable boson self-energy — in*(£i, ta), and the Auger 
component of the ground state boson self-energy —iU g (ti,t2). 
In the upper most diagram the + sign in the indices means 
that there are two separate diagrams, one with the k and * in- 
dices and one with the q and g indices, both of which con- 
tribute in an additive way and have thus to be summed. Fur- 
thermore in the diagram below the upper most one the sub- 
scripts k and hold for 11* and the subscripts q and 1 hold 
for II g . In all diagrams a double line indicates a full propaga- 
tor whereas a single line denotes an unperturbed propagator. 



with 



h 2 ^ fc /<r 

k/q 



o/i 



kq 



(22a) 



(22b) 



xGf\t 2 ,h)Gf (ti.ta), 



and G^£ and B^J X denoting, respectively, the contour- 
ordered Green functions for a free/ valence band electron 
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and an auxiliary boson. 



In the self-energies (21) the two reaction channels ([I]) 
and ^ are separated. Every term involving or 
refers to the RCT channel and every term containing ap- 
pertains to the Auger channel. Due to the dressed Green 
functions the two channels are however coupled. The 
Green functions of the auxiliary bosons contained in op 
and Oq ensure that the two tunneling processes contained 
in ([I]) are resonant. Physically, the auxiliary bosons simu- 
late the action of intra-molecular correlations which kick- 
in when an electron hops to-and-fro the molecule. 



Using the Langreth-Wilkins rules for analytic contin- 
uation (see Refs. I3T1 and 1551 and Appendix [A} we obtain 
from the self-energies (21 1 the set of Dyson equations 



given in Appendix[B] The components of these equations 
arising from the RCT terms are equivalent to the ones in 
Ref. I3H but with two bosonic pseudo-particles instead of 
one and an energy shift caused by the auxiliary bosons. 
The set of Dyson equations (Bl| contains terms which 



violate the constraint (14). Before physically meaningful 



information can be extracted the Dyson equations have 
to be projected onto the physical subspace defined by the 
constraint (14). 



The procedure to achieve this is originally due to 
Langreth and Nordlander and has been outlined several 
timeS^HH] it is based on inspecting the order of the 
Green functions in the conserved pseudo-particle number 
Q. The retarded functions G_ and B R , g are proportional 
to Q° while the lesser Green functions and B^, g are 

proportional to Q 1 . Thus, we have to omit any terms of 
higher order than Q° from the retarded self-energies and 
any terms of higher order than Q 1 from the lesser self- 
energies. This approach is not an additional approxima- 
tion but an exact projection enforced by the constraint. 



Before carrying out the projection we split off the 
Green functions' oscillating factors by means of the de- 



compositions^ 

Gt /R/A (t, = G< /R/A (t, f) e-*X5 dtl 



(23a) 



B 



</R/A 



t\ _ u<l R / A 



(t, if) = B 



(M')e-5U' dtl£ * (tl) ) (23b) 



B</ R / A {t,t') = j?</«M(t,i')e-*# d * ie «< tl > , (23c) 
and 

G R (t, if) =-iQ(t- t') g_ (t, if) , (24a) 

G A (t, if) = i B(lf - t) g_ (t, if) , (24b) 

B R /g (t, t') = -i Q(t - t') K /g (t, if) , (24c) 

B A /g (t,t')=ie(t' -t)K /g (t,t') . (24d) 

Using the definition of the retarded and advanced Green 
functions we find the following relations^ 

g-(t,t) = b* /g (t,t)=l, (25a) 

g-(t,t')= [g-(t',t)]*, (25b) 

K /g (t,t') = [K /g (t',t)]* . (25c) 

Within the Dyson equations the oscillating terms emerg- 



ing from Eqs. (23) will be absorbed in the functions ct5, 
a- and which are defined by 

1 kq J 

a|(t,0 = ^(* I Oe"* / *' d ' lle - (tl) " e * (il)] . (26a) 
of (M') =Z>(t,t')e-iti dH t e -^- e ^ , (26b) 
a>(t,t') =a?(t,i')e"* J ' t ' dtl[e * (tl) " e9(il)] . (26c) 

kq y ' kq K ' v ' 

The terms a2", a$, and vanish identically due to 

the initial conditions n (£o) = n g(^o) = 1 an< i "-i(io) = 
n q(to) = since 

ff|(t,0~ (l-»j(*o))(l-no(*o))=0, (27a) 

(27b) 
(27c) 



trf(t,lf) ~ n f (t )m(t ) =0, 



Employing the Langreth-Nordlander projection together 
with the relations p3|, pi|, J26| and p7l) the set of 



Dyson equations (Bl) takes the following form 



|b<(«,«') 
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dt 1 S : |(t,t 1 )6 ff (t,ti)Gf(ti,t / )+ / dti?|(t,ti)B<(t,ti) ff _(t 1 ,t / ) 



dti 



ct| (t ! , t) 5 _ (t, tj) B< (*!,£') + i5^(t, i x ) 6 S (i , t x ) B< (h , t') 



dti 



5> t) G< (t, ti) b g (ti, t') + ioF> f (t a , t) B< (t, t! ) fe g (tl , t') 



Q(t-t')~g_(t,t') = - I dtiS£(t,f 1 )& ff (i,ti)ff_(ti,0 , 
e(t-t')J^*(M') = I d/, 

e(t-t')^MM') = o 



(28a) 
(28b) 
(28c) 
(28d) 
(28e) 
(28f) 



The time evolution of the molecular occupation num- 
bers ri- , n* and n g can be calculated from the following 
relations^ 



(29a) 
(29b) 
(29c) 



dn_(t) 


dG<(t,t') 


dG<(t,t') 




dt 


dt 


t=v dt' 




dn* (t) 


9B<(t,t') 






dt 


dt 




t=t 


dn ff (t) 


dB<(t,f) 


, dB<(t,t') 
H 2 




dt 


dt 




t=t 



Equations (28) and (30 1 constitute the final projected 
set of Dyson equations that determines the dynamics 
of the system within the subspace Q = 1. The rate- 
cquation-like structure of these equations is already 
evident. The Dyson equation for the lesser Green 
function of the negative ion, Eq. (30a I, for instance, 



contains a production term proportional to ov and B< 

and a loss term proportional to a$ and Gf . These terms 
relate to the production and loss of negative ions by the 
RCT electron capture and release reaction, respectively. 



Hence, we also need the adjoint Dyson equations of the 
lesser functions which can be calculated in the same man- 
ner. The result reads 



dt 



7 G<(t,t') 



d_ 

dt: 



7 B<(t,t') 



t' 



dt, G^t^ih^bgihJ) 



Ah 5 _(t,ti)a|(ti,t').B<(ti,t') , 



(30a) 



d*x B<(tM)^{t'M)g-{ti,-t') 



dti S<(t,ti)i??_,(ti,t')& g (ti,t') , 



(30b) 



SEMI-CLASSICAL APPROXIMATION 



The strongly oscillating factors of the projected set of 
Dyson equations are contained in the functions aS, <r|', 
and <Tg_. If these functions are strongly peaked along the 

time diagonal t = t' , we can apply a saddle point approx- 
imation to the integrals in ( 28 ) and ( 30 1 . For instance, 



d^ af(t,ti)B<(t,h)g-(ti,i/) 



B<(t,t)g-(t,t / ) I dtia|(t,ti) 

-OO 



(31) 



d_ 
dt- 



dti &»(*>*i)3f (i'.fOG^ti.t') 



Ahbgit^i^St'^B^hX) 



(30c) 



The validity of this approximation, which is also known 
as the semi-classical approximatiorPS will be demon- 
strated in Sec. IVII 

Within the saddle-approximation the projected Dyson 
equations for the lesser Green functions become 
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d 

dt" 
d 



G<(t,t')n-b g (t,t)G<(t,t') / dfic£(t,fi) + B<(f,t) fl _(t,t') / &tia<{t,tx) 



t' 



^G<(t,t') ^ -G<(t,t')b g (t',t') I dti*$(ti,ir) + g-(t,1?)B<tf,l!) I dt x u<{hj) 

d 



dt 

d_ 

dt 

d 



B<{t,t r )^-g_(t,t)B<{t 7 t r ) J dti ff^Cti,*) - & fl (t,t)S<(t,0 / dhia^h) 
i 1 

7 5<(t,i') » -Bf{t,t')g-{t',t') f dti - B<(t,t/)b g (t',t'\ / dti i%(h,t' • 



dt 9 

d_ 
dt/ 



B<{t,t')^G<{t,t)b g (t,t r ) I dt 1 a^{t 1 ,t)+B<{t,t)b g (t,t') I dh ia> q ih,t) , 



— OO 



7 B<(t,i/)^b g (t,i/)G<(t',t') I dilo^>(^ , ,tl) + 6 fl (t 1 i / )s<(f / > ^ , ) / dMSf^i) ■ 



(32a) 
(32b) 
(32c) 

(32d) 

(32e) 
(32f) 



Using Eqs. ( 29 ) we then arrive at a set of rate equa- expressions for these two rates - as obtained from the 
tions for the occupancies of the molecular pseudo-particle semi-empirical model - can thus be looked up in our pre- 
states, 



vious worw^ 21 . The Auger rate Ta introduced in (34c 



dn-(t) 

dt 
dn* (t) 

dt 
dn g (t) 

dt 



-ri(t)n_(t) + r (t) n»(t) 

-r (t) n*(t)-r A (t) «•(«) 
ri(t)n_(t) + r A (t)»i*(t) , 



where the rates are given by 

r (i) = / d^ 2K{c£(t,ti)} , 

«/ — OO 

r x (t) = / dt 1 23?{a>(t,t 1 )} + r n , 

«/ — OO 

I\i(*) = T dt 1 2K{»oF>(t,t 1 )} . 

— OO 



(33a) 
(33b) 
(33c) 

(34a) 
(34b) 
(34c) 



however has not been calculated before. Within the semi- 
empirical model it is given by 



d<? r 25R 
dq# 



fi 2 (Aq) 3 (A/fc) 3 

dgv 



dti 



dk r 



2tt 



dks / dk v (35) 



x sin(g^) kf, sin(A^)rg f (*i) 



with 



r Sf (t) = V^ ff (t)eK/o dt iN(tO+M*i)-^(ti)-^] . (36) 



The multi-dimensional integral in ( 35 ) can be calculated 
efficiently using the approximations outlined in Ref. I27I 
We will now seek an analytic solution to the coupled 
the negative ion by adding the natural decay rate T n — rate equations (33). As a starting point we first take a 



Note, in Eq. (34b) we incorporated the natural decay of 



1 jr n on the right-hand side. 

Similar to what Langreth and coworkers did in the con- 



step back and consider the isolated decay channels of res- 
onant electron capture, resonant electron emission, and 



text of the neutralization of atomic ions^S we have thus Auger de-excitation. Singling out the individual reac- 
reduced a complicated set of Dyson equations - Eqs. ((28|) tions in (33) we obtain 

dn^\t) 



and ([30]) describing the de-excitation of a metastable 
molecule via the simultaneous action of the RCT chan- 
nel and the Auger channel ([2| - to an easy to handle 
system of rate equations ( 33 1 . The reaction rates ( 34 ) en- 
tering the rate equations are linked to quantum-kinetic 
quantities and thus related to the semi-empirical model. 

The rates r and defined in (34a) and (34b) are 
equal to the rates employed in Ref. 28. Mathematical 



dt 



(*) 



dt 

dn^jt) 
dt 



= -r (t) n! 0) (t) 



-r x (t) n W (i) 



-T A (t)ni A \t) . 



(37a) 
(37b) 
(37c) 
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,(o) 



and 



The superscripts (0), (1), and (A) identify the isolated 
resonant electron capture, resonant electron emission, 
and Auger de-excitation, respectively. Since the channels 
are isolated, each of the decay equations (37) comes with 



completely determined by the occupancies n 

n{ A \ Moreover, when the Auger channel is disabled by 
setting T A (t) ee 0, Eqs. Q, @ and (El reduce to the 



an analogous equation for the species that is produced. 
For instance, accompanying (37a) is the equation 



dn W (t) 
dt 



(38) 



The time derivatives of nl°' and differ however only 



in sign. Hence, n_ is given through the conserva- 
tion of particles as = 1 — (valid when the chan- 
nels are isolated). Consequently, the additional equa- 
tions of type ( 38 ) do not contain any additional infor- 



mation and can be omitted. Using the initial condi- 



tion nf^ (to) — nf 1 ' (to) — 1 the system ( 37 ) can be solved 
straightforwardly. The result is 



rate equations derived by intuitive means for the isolated 
RCT channel. Hence, the quantum-kinetic treatment 
justifies a posteriori the intuitive approach taken by us 
in Ref. M 

We now turn to the spectrum of the emitted electron. 
While the evolution of the q states has not been con- 
sidered explicitly in our quantum-kinetic calculation, the 
occupancy of these states can nevertheless be extracted 
from the solution of (33). 

From the reactions (|1| and ^ it is obvious that the 
probability for emitting an electron n e (t) is equal to the 
occupancy of the ground state n g (t) as every ground 
state molecule must have resulted from the reaction chain 
and, hence, must be accompanied by an emitted elec- 
tron. Consequently, the evolution of n e (t) is governed by 
Eq. (33c). Due to the image potential, however, not ev- 



n { °\t) =e 



n ( l\t)=n ( }\t>)e-ti^(t) 



ni A) (t) = e 



(39a) 
(39b) 
(39c) 



Now we are in position to use these occupancies to 
calculate the solution of the full, coupled system of rate 



equations ( 33 1 . First we consider the equation for 



Using the initial condition n*(io) = 1, Eq. (33b) can be 
solved by separation of variables and yields 

nS) = e -4<«x(ra(t0+ri<t0) = n (°) {t) n W {t ) . (40 ) 



ery emitted electron can escape the surface. In particular, 
the escape is only possible when the emitted electron's 
perpendicular kinetic energy = e?" cos 2 (q$) is higher 
than the absolute value of the image potential Vi at the 
position of emission. The latter can be approximated by 
the position of the molecule's center of mass at the time 
of emission. 

To incorporate the image potential effect, we adopt a 
two step strategy. As a start we introduce the spectral 
rates Q^s"^ ,t) and QA(^,t) which are not restricted by 
the image potential effect by writing 



T 1/A (t)= / defg 1/A (ef,t) 
Jo 



(44) 



To solve (33a) for the occupancy of the negative and afterwards we let Qi/a — > Qi/a with 



ion state we first multiply this equation by a factor 
exp(f t ' d^Ti^)) and rearrange the terms to obtain 



dt 



■(*) 



ft dt 2 ri(t 2 )\ ^ ,,n /.' dt 2 r 1 (t 2 ) 



(41) 

Relabeling then t as t\ and integrating the equation 
from t\ — to to t\ — t while minding the initial condition 
ri_(<o) = yields after a further rearrangement 



Qi/ A W,t) 

x e(Vi(z R (t)) 



dq# / dg ¥ 

n JO 

x , A 2 Qi/ A (ef,t) 



L-9 



(45) 



dq-gdq, 



An explicit expression for the spectral RCT emission 
rate Q\ has been given in Ref. [28l The spectral Auger 



_(*)=/ dta T (tOn, (h ) e- A dt2Fl(t2) 



= / dtj 



dn^jh) 
dii 



rate qa may be calculated from Eq. ( 35 ) by stripping out 
the q r integral and multiplying the result by m e /h 2 q r . 
Introducing the spectral decomposition of the rates in 

with n e 



(42) Eq. ( 33c ) and identifying n s 



we obtain 



Finally, the occupancy of the molecular ground state 



dn(t) 
dt 



n g , that is, the solution of Eq. (33c), is given through the 



def g^ef^n-it) 



+ defQ A (e?,t)n*(t), 



particle conservation property of the full system ( 33 1 , 

n g (t) = l-n*(i) -n_(i) . (43) 

Note, the molecular occupancies satisfying the com- 
bined rate equation scheme Eqs. (40), (42), and (43) are 



(46) 



where n e denotes the probability for emitting an electron 
that can escape from the surface. Integrating over the 
time argument with the initial condition n e (to) = and 
taking the derivative with respect to e°? we find for the 
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spectrum of the emitted electron at time t 
dn P 

t 



def 



dti gi(e§P ,ii)n_(*i) 



(47) 



to 



The secondary electron emission coemcient 7 e , that is, 
the probability for having emitted an electron after the 
collision is completed, is given by 



7e 



n e (oo) 



(48) 



The occupancies of the molecular pseudo-particle 
states ( 40 1 , ( 42 ) , ( 43 ) and the spectrum of the emitted 



electron ( 47 1 are the main result of this work. The occu- 



pancies fully characterize the temporal evolution of the 
de-excitation of a metastable N2( 3 E+) molecule at a sur- 
face when both the RCT and the Auger channel are open. 
The ingredients required as an input, the occupancies 
arising from the isolated processes , n_ , ni A ^ and 
the image potential adjusted spectral rates Q\, qa, can 
be obtained from the quantum-kinetic calculation and 
thus from the semi-empirical model. Since the param- 
eters of the model are assumed to be fixed, either by 
experiment or by quantum-chemical calculations, there 
is no free parameter in our theory. 



VI. NUMERICAL RESULTS 

In this section we present numerical results based on 
the semi-classical equations of Sec. |Vj We consider the 
particular case of a diamond surface and restrict our in- 
vestigations to normal incident with a molecular kinetic 
energy of 50 meV. The turning point of the molecule's 
traj ectory is then 4.4 Bohr radiPn As in our previous 
workpZlSS] we treat only the two principal orientations of 
the metastable N2( 3 E+) molecule: molecular axis per- 
pendicular to the surface and molecule axis parallel to 
the surface. Furthermore, we omit the surface induced 
decay channel by setting Ti(t) = T n . As our previous in- 
vestigations showed, this is an excellent approximation^. 

The numerics necessary to calculate within the semi- 
empirical model the Auger matrix element V^- has been 
described in Ref. [27J Utilizing the fact that for the low 
collision energies we are interested in the turning point 
of the molecule is far outside the solid as well as the fact 
that the LCAO wave functions are strongly localized on 
the molecule whereas the wave functions of the solid and 
free electron are bounded in the mathematical sense we 
split-off the time-dependence of the matrix element and 
integrate the rest by an interpolative grid-based Monte 
Carlo scheme. The tunneling matrix elements Vt and 
can be calculated within the semi-empirical model partly 
analytically and partly numerically. The required tools 
are given in Ref. |2"51 



(b 
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FIG. 4. Variation of the real (solid line) and imaginary 
(dashed line) part of a^(ti,t2) (upper panel) and a%Jti,t2) 
(lower panel) as calculated from ( |22[ ) along the anti- 
diagonal ti = — 12 = t. The molecule's axis was aligned per- 
pendicular to the surface and the molecule's kinetic energy 
was fixed to 50meV. The behavior for negative time argu- 
ments t is omitted since the real (imaginary) part of both 
functions is symmetric (anti-symmetric) with respect to the 
time diagonal. Note that the time t is dimensionless. It re- 
lates to the physical time t phyB via t p h ys = a B At/2vo- 



First, we investigate the validity of the semi-classical 
approximation. As outlined above, approximations of the 
form (31 1 are only acceptable if the functions crg^i^) 
and a^Jtx^ti) are sufficiently peaked on the time diag- 
onal t\ =t%. In order to validate this assumption for 
our case, we plot in Fig. [4] for the perpendicular orienta- 
tion the variation of these two functions along the anti- 
diagonal t\ = — *2- The plots can be generated directly 
from the definitions ( |22[ ) and represent profiles with re- 
spect to the time diagonal. 

For both functions c?g(ti,t2) an d £2) the real 

part has its maximum on the time diagonal whereas the 
imaginary part vanishes on the diagonal itself but ex- 
hibits the highest value very close to it. When the time 
separation from the diagonal is enlarged, both functions 
decrease in an oscillating way. For 3v the amplitude de- 
creases to about 10% over a dimensionless time inter- 
val of At m 0.05. This relates to a physical time span 
of Atphys = a,BAt/2v 2.25 fs and to the motion of 
the molecule along a distance of 0.025 as- For cfg^. the 
fall-off is even more drastic. The behavior for shifted 
anti-diagonals as well as for parallel orientation is very 
similar and hence not shown here. Altogether, we can 
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incoming 



outgoing 




FIG. 5. Time dependence of the occupancy of the ground 
state molecule n g , the metastable molecule n*, and the neg- 
ative ion n_ in the parallel (solid line) and perpendicular 
(dashed line) orientation. The molecule's kinetic energy was 
fixed to 50 meV. The dotted line represents half filling of the 
respective level. The y-axis is logarithmic and the time is de- 
noted in the dimensionless units of Fig. [4] The curves were 
calculated from Eqs. (401, (42 1 and (|43[|, respectively. The in- 



coming and outgoing branches of the trajectory are indicated 
at the top of the diagram. 



conclude that with respect to the macroscopic motion of 
the molecule the functions trj and are indeed suf- 
ficiently peaked on the time diagonal. Thus, the semi- 
classical approximation is valid in our case. 

We now turn to the occupancies of the molecular 
pseudo-particle states, that is, the probabilities with 
which the molecular configurations involved in the de- 
excitation process appear in the course of the scat- 
tering event. The time dependent occupancy of the 
ground state n g , the metastable state n*, and nega- 



tive ion state n_ can be calculated from Eqs. (40), (42) 



and ( 43 ) , respectively. The results are depicted in a semi- 



logarithmic plot in Fig. [5j Inspection of the curves re- 
veals that even close to the surface the occupancy of the 
negative ion state is rather low. Hence, the metastable 
molecule is almost immediately converted into a ground 
state molecule. In Fig. [5] this fact is recognizable at the 
crossing point of the n* and n g curves which occurs at 
approximately half filling of both levels. The low occu- 
pancy of the negative ion state is caused by the high 
efficiency of the natural decay channeP^ and not by the 
Auger channel destroying the metastable molecule, which 
is the generating species of the negative ion. In fact, it 
is the other way around and in order to substantiate this 
claim, we investigate next the relative efficiency of the 
RCT and Auger channel by considering the respective 
reaction rates. 

Figure [6] shows the rates of resonant electron cap- 
ture r and Auger de-excitation Ta- For both chan- 
nels the rates are highest at the molecule's turning point 
(approximately 4.4 as) which is the point of small- 




FIG. 6. Variation of the rates of resonant electron capture To 
(upper panel) and Auger de-excitation Fa (lower panel) for 
the parallel (solid line) and perpendicular (dashed line) ori- 
entation at a kinetic energy of 50meV. The time is denoted 
in the dimensionless units of Fig. [4] 



est molecule-surface separation and strongest molecule- 
surface interaction. When the molecule-surface distance 
is increased, the rates decrease exponentially. The RCT 
channel's rate is about two orders of magnitude higher 
than the Auger channel's rate. Consequently, the RCT 
channel captures surface electrons much more efficiently 
than the Auger channel. In fact, the RCT channel 
is so effective in capturing electrons that it under-runs 
the Auger channel by destroying its starting basis, the 
metastable state. As a result, in the combined two- 
channel system the Auger channel's performance is sig- 
nificantly diminished as compared to the isolated Auger 
reaction. 

This conclusion may be verified by considering the 
term in the rate equations which is responsible for the 
production of the ground state molecule by an Auger de- 
excitation. It is given by (see Eqs. (|33c| and (40)) 



T A (t) n*(t)=T A (t) nl A \t)n<?\t) 



(49) 



Here the factor n^\t) is only present in the combined 
two-channel system but not in the isolated Auger system. 
Without explicit proof but based on numerical observa- 
tions we note that the term (t) is almost identical to 
the combined occupation n*(i) depicted in Fig. [5} Hence, 
in the combined system the Auger channel's ground state 
production term ( 49 ) is strongly suppressed already in 



the incoming branch of the trajectory. 
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FIG. 7. Energy spectrum of the emitted electron in paral- 
lel (upper panel) and perpendicular (lower panel) orientation 
calculated from ( |47[ l at t = oo for a kinetic energy of 50meV. 
In both panels the dotted line specifies the isolated RCT spec- 
trum (obtained by setting qa = 0), the dashed line denotes 
the isolated Auger spectrum (obtained by setting Qi =0), and 
the solid line represents the combined two-channel spectrum. 



Finally, we turn to the energy spectrum of the emit- 
ted electron. Figure [7] depicts the emission spectrum at 
t = oo for the combined two-channel reaction as well as 
for the isolated reaction channels. The latter can be ob- 
tained by setting in ( [47] ) qa = or g± = 0. As can be 
seen, the isolated RCT spectra exhibit a strong peak at 
about 1.5 eV and slowly drop off for higher energies. The 
isolated Auger spectra, on the other hand, monotonously 
increase until approximately 2.8 eV and then immedi- 
ately fall off. The low energy cut-off of all curves is 
due to the trapping of the emitted electron in the im- 
age potential close to the surface when its perpendicular 
energy is too low. The combined spectra almost equal 
the respective isolated RCT spectra. Only in the range 
from 1.5 eV to 2.5 eV are the combined spectra slightly 
increased with respect to the RCT curves. This minor 
enlargement is due to the Auger channel and supports 
our previous finding that the RCT channel dominates 
the Auger channel. 

The combined spectra in Fig. [7] are different from the 
simple addition of the isolated spectra. This behavior is 
caused by the unified treatment of the RCT and Auger 
reaction channels. The effect would be even more pro- 
nounced for molecular species forming stable negative 
ions. Here the resonant electron emission would be al- 
most completely blocked as the surface induced decay is 



TABLE I. Secondary electron emission coefficients in parallel 
(7I) and perpendicular (7^) orientation at a kinetic energy 
of 50 meV. 



always very weak. The resonant electron capture, how- 
ever, would be still very efficient in destroying the initial 
species. Consequently, the spectrum of the emitted elec- 
tron would resemble the Auger spectrum in shape but 
would be strongly decreased in magnitude. 

The secondary electron emission coefficients are given 
by the area beneath the curves in Fig. [7] and are sum- 
marized in Table [I] In accordance with our previous ob- 
servations the emission coefficients are not changed sig- 
nificantly by the inclusi on of the Auger channel. A simi- 
lar result was found by Stracke et al.'^for N 2 ( 3 £+) de- 
exciting at a tungsten surface. Their experimental mea- 
surements imply that only about 10% of the secondary 
electron emission coefficient is made up by the Auger 
channel. 



VII. CONCLUSIONS 

We constructed in this work a semi-empirical gen- 
eralized Anderson-Newns model for secondary electron 
emission due to de-excitation of metastable N2( 3 £+) 
molecules at dielectric surfaces. The model treats Auger 
de-excitation and the two-step resonant charge transfer 
process, where the Nj( 2 II fl ) ion acts as a relay state, 
on an equal footing. It reduces the molecular projec- 
tile to a two-level system representing the molecular or- 
bitals which change their occupancies during the reac- 
tion and treats the surface as a simple step potential 
confining the electrons of the solid. By construction, 
the semi-empirical model is not restricted to a partic- 
ular projectile-target combination. Having applications 
of the model to charge-transferring processes at plasma 
walls in mind, where a great variety of different projectile- 
target combinations occurs, we consider this as a real ad- 
vantage. Another advantage is that the semi-empirical 
model separates the many-body theoretical description 
of the non-interacting projectile and target from the 
quantum-kinetic treatment of the scattering process. The 
former is simply encapsulated in the parameters of the 
model Hamiltonian and the latter is performed by Green 
functions. This is particularly advantageous in cases 
where the surface scattering event is studied primarily 
because of its connection to the physics of quantum- 
impurities. 

For the semi-empirical model to work a method was re- 
quired to assign and control the energies of the two-level 
system in accordance to the reaction channels, that is, to 
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have the two-level system describing all three molecular 
configurations involved in the de-excitation process: the 
metastable N 2 ( 3 £+) molecule, the negative ion N^~( 2 II ff ), 
and the molecular ground state N2( X S+). We showed 
how this can be done with projection operators and aux- 
iliary bosons. As a result, both the resonant tunnel- 
ing and the Auger channel could be cast into a sin- 
gle model Hamiltonian which, with the help of pseudo- 
particle operators, could then be made amenable to a 
diagrammatic quantum-kinetic calculation. Using the 
self-consistent non-crossing approximation for the self- 
energies and a saddle-point approximation for the time 
integrals in the self-energies we finally derived from the 
Dyson equations for the propagators of the molecular 
pseudo-particles a set of rate equations for the probabili- 
ties with which the molecular configurations contributing 
to the de-excitation process can be found in the course 
of the scattering event. Without the Auger channel, the 
system of rate equations reduces to the one postulated 
by us before on intuitive grounds for the RCT channel 
alone^. The present work justifies therefore this reason- 
ing a posteriori. 

For the particular case of a diamond surface we veri- 
fied the validity of the semi-classical approximation and 
investigated for a collision energy of 50meV the inter- 
play of the resonant tunneling and the Auger channel. 
In particular, we analyzed the temporal evolution of the 
probabilities with which the projectile is to be found in 
the N 2 ( 3 S+), the N^( 2 II g ), or the N 2 ( 1 S+) state and 
explicitly calculated the rates for electron capture due 
respectively to tunneling and Auger de-excitation. We 
also obtained the spectrum of the emitted electron and 
the secondary electron emission coefficient 7 which are 
the two quantities of main importance for the modeling 
of gas discharges. Our results indicate for a diamond sur- 
face and a kinetic energy of 50 meV the resonant tunnel- 
ing channel clearly dominating the Auger channel. The 
contribution of the Auger channel to the secondary elec- 
tron emission coefficient lies only in the range of a few 
percent. The overall 7 coefficient is on the order of 10 _1 
in agreement with what has to be typically assumed to 
make kinetic simulations of dielectric barrier discharges 
reproduce the properties of the discharge. 

With minor modifications the semi-empirical model 
and its quantum-kinetic handling leading to the easy 
to use set of rate equations can be adopted to other 
plasma-relevant charge-transferring surface collisions as 
well. Thereby we hope to replace at least partly the 
rules of the thumb which are often needed to characterize 
secondary electron emission due to neutral and charged 
heavy plasma species hitting the plasma wall. 
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Appendix A: Langreth-Wilkins rules 

The Langreth-Wilkins rules^ are a powerful tool for 
the analytic continuation of propagators defined on a 
complex time contour onto the real time axis. Their ex- 
plicit form depends on the initial definition of the Green 
functions. Unfortunately, there is no common agreement 
about the usage of i factors. Moreover, rules published in 
the past sometimes contained typographic errors^. Due 
to these reasons we list the explicit form of the Langreth- 
Wilkins rules used in this work. The rules can be derived 
in the standard wa^ * 35 * 40 * using however the definitions 



( 17 )-( 19 ) for the Green functions. In the following F and 



B denote fcrmion and boson propagators, respectively. 

To analytically continue the boson-like fermion- 
antifermion pair 



B(t,t') = F l (t,t')F 2 (t',t) 



we utilize 



B>(t,t>) 
B<(t,t') 

B R (t,t') 



B A (t,t>) = 



= iF>(t,t')F<(t',t) , 
= iFf(t,t')F>(t',t) , 

F<(t,t')F A (t',t) 
+ F R {t,t')F<{t',t) 
F<( tl t')F*(t\t) 
+ F A (t,t')F<(t',t) 



For the fermion-like fermion-boson pair 

F{t,t') = F 1 {t,t')B 1 {t,t') , 

the following rules hold 

-iF>(t,t')B>(t,t') , 
-iF<(t,t')B<(t,t') , 

-i[F R (t,t')B<(t,if) 
+ F>(t,t')B R {t,t') 
F A (t,t') = - i [F 1 A (t,t')B<(t,t') 

+ F>(t,t / )B A (t,t') 



F>(t,t') 
F<(t,t') 

F R (t,t') 



(Al) 



(A2a) 
(A2b) 

(A2c) 



(A2d) 



(A3) 



(A4a) 
(A4b) 

(A4c) 



(A4d) 



Johannes Marbach was funded by the federal state of 
Mecklenburg- Western Pomerania through a postgradu- 
ate scholarship. In addition this work was supported by 



The boson-like boson-boson pair 

B{t,t')=B 1 {tX)B 2 {t,t r ) , 



(A5) 
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can be analytically continued by 



B>(t,t') 
B<(t,t') 

B R {t,t') 
B A (t,t>) 



iB>(t,t')B>(t,t') , 
iB?(t,t')B<(t,t') , 

B R (t,t')B>(t,t') 
+ B<(t,t / )B R (t,t / ) 
B A (t,t')B>(t,t') 
+ B<(t,t')B A (t,t') 



(A6a) 
(A6b) 

(A6c) 



(A6d) 



Finally, to analytically continue the boson-like boson- 
antiboson pair 



B{t,t')=B x {t,t')B 2 {t',t) , 



we use 



B>(t,t') 
B<(t,t') 

B R (t,t') 
B A (t,t>) 



-iB>{t,t')B<{t',t) , 
-iB<(t,t')B>(t',t) , 

-i[B<(t,t')B A (t',t) 
+ B R (t,t')B<(t\t) 

-i[B<{t,t')B R {t',t) 
+ B A (t,t')B<(t',t) 



(A7) 



(A8a) 
(A8b) 

(A8c) 



(A8d) 



Furthermore, for the analytic continuation of the con- 
tour integrals within the Dyson equations we also need 
to project terms of the form 

D{t,t') = J^dh £>i(Mi) D 2 (t u t') , (A9) 

where D, D\ and D 2 are either fermion-like or boson-like. 
This can be accomplished by the rules 



£><(*, i') 



D>(t,t') 



dti 



D R {tM)D<{h,t') 
-Dfit^DiihJ) 



D R {tM)D>(t u t') 
+ D>(t,t 1 )D A (t 1 ,t') 

/OO 
dh D R (t,h) D?(ii,f) , 
-OO 

/OO 
dh D A {tM)D A {t 1: t') ■ 
-OO 

Appendix B: Dyson equations 



(AlOa) 

(AlOb) 

(AlOc) 
(AlOd) 



In this appendix we summarize the Dyson equations 
for the analytic pieces of the molecular Green functions 
G_, and B g as obtained by an application of the 
Langreth-Wilkins rules of appendix|A"|to the self-energies 
shown in Fig. [3] 

The lesser Green functions satisfy 



16 



'.d_ _ e_(t) 
l dt h 



G<(t,t') = / dti c7>(t,t 1 )Sf(t,t 1 )G<(t l! + ^(*,*i)S<(t,t 1 )G<(ti,t / ) 



+ a> (t, t x ) Bf (t, t x ) Gf (t x , t') + <r5(t, t x ) B< (t, t x ) G<(h,t') 



(Bla) 



dti 



a5 (t, tl ) B< (t, h) (h ,t')+ <r< (t, ti) B< (t, h) G± (h , t') 



.d_ _ ejfy 

l dt h 



B< (t,t')= dti k< (t i , t) G R \{t , ti ) B< (ti , t') + <r^(ti , t) Gf (t, ti) fl< (ti , t') 

J —oo 

+ ia>Jfi, ti) Sf (t, ti) B<(ti, t') + iof^t, ti) S<(t, ti) B<(ti,t') 



(Bib) 



St 



dti 

J — OO 

Bf{t,t')= f dti 
+ / dti 

while the retarded Green functions obey 
d e_(t) 



(ti, t) G< (t, ti) B^ 1 (ti, t') + » CT < (t, ti) B< (t, ti) B, A (ti, t') 



kq 

a< (ti, t) G« (t, ti) B< (ti , t') + tr^(ti , t) G< (t, ti) S< (t x , t') 
+ ^| f (ti,t) Sf (t, ti) B<(ti, t') + ia^(ti,t) S<(t, ti) B<(ti,t') 



(Blc) 



a} (ti , t) G< (t , t x ) B g A (t i , t') + ia> (h ,t)B<(t,h) (t i , t') 



at 



G R {t 1 t')=5(t-t') + / dt 1 [a>(t,t 1 )B R (t,t 1 )G R (t 1 ,t / )+af(t 1 t 1 )B<(t,t 1 )G R (h,t') 

+ g> 3 (t, ti) B?(t, ti) G R {h,t') + «r5(t, ti) B< (t, ti) G«(ti, t') 



(Bid) 



._a _ e*(t) 

at ~1T~ 



. d e g (t) 
dt h 



B R {t,t>)=8{t-t')+ I dti c7<(ti,t)G^(t,ti)Sf(ti,t')+^(ti,t)G<(t,ti)Sf(ti,t') 



ti) S« (t, ti) Bf (ti, t') + ia£_.(t, ti) B<(t, ti) Bf (ti, t')] , 

(Ble) 



B*(t,t')=5(t-t')+ I dti af{t ll t)G R {tM)Bf{t ll t') + a^{t ll t)G<(tM)Bf{t ll t') 

+ ia<Jt u t) B R (t, h) S«(ti,t') + ^(ti,t) B<(t, ti) B R (ti, t') 



(Blf) 



The greater and advanced Green function can be obtained from the definitions ( 19 ) and the symmetry relations ([20] 
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